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^ ■ Abstract. We present a map of 511 keV electron-positron annihilation emission, based on data accumulated with 

\Q ' the SPI spectrometer aboard ESA's INTEGRAL gamma-ray observatory, that covers approximately ~ 95% of 

the celestial sphere. Within the exposed sky area, 511 keV line emission is significantly detected towards the 
galactic bulge region and, at a very low level, from the galactic disk. The bulge emission is highly symmetric and 
is centred on the galactic centre with an extension of ~ 8° (FWHM). The emission is equally well described by 
models that represent the stellar bulge or halo populations. The detection significance of the bulge emission is 
~ 50cr, that of the galactic disk is ~ 4<7. The disk morphology is only weakly constrained by the present data, 
being compatible with both the distribution of young and old stellar populations. The 511 keV line flux from the 
bulge and disk components is (1.05 ± 0.06) x 10 -3 ph cm _2 s _1 and (0.7 ± 0.4) x 10 -3 ph cm _2 s _1 , respectively, 
Oh' corresponding to a bulge-to-disk flux ratio in the range 1 — 3. Assuming a positronium fraction of f p = 0.93 this 

translates into annihilation rates of (1.5 ± 0.1) x 10 43 s _1 and (0.3 ± 0.2) x 10 43 s _1 , respectively. The ratio of the 
bulge luminosity to that of the disk is in the range 3 — 9. We find no evidence for a point-like source in addition to 
the diffuse emission, down to a typical flux limit of ~ 10~ 4 ph cm -2 s -1 . We also find no evidence for the positive 
latitude enhancement that has been reported from OSSE measurements; our 3cr upper flux limit for this feature 
is 1.5 x 10" 4 ph cm 2 s . The disk emission can be attributed to the /3 + -decay of the radioactive species 26 A1 
and 44 Ti. The bulge emission arises from a different source which has only a weak or no disk component. We 
suggest that Type la supernovae and/or low-mass X-ray binaries are the prime candidates for the source of the 
galactic bulge positrons. Light dark matter annihilation could also explain the observed 511 keV bulge emission 
characteristics. 
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1. Introduction 

Since the first detection (Johnson & Haymes 1973) and 
the subsequent firm identification (Leventhal et al. I1978() 
of the galactic 511 keV annihilation line, the origin of 
galactic positrons has been a lively topic of scientific 
debate. Among the proposed candidates for sources of 
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positrons figure cosmic-ray interactions with the inter- 
stellar medium (Ramaty et al. I197U|I . pulsars (Sturrock 
1971), compact objects housing either neutron stars or 
black holes (Ramaty & Lingenfelter 1979), gamma-ray 
bursts (Lingenfelter & Hueter H984II . (light) dark matter 
(Rudaz & Stecker ITTjggl Boehm et al. l2UUljl . and stars 
expelling radioactive nuclei produced by nucleosynthesis, 
such as supernovae (Clayton 1973), hypernovae (Casse 
et al. 12004( 1 . novae (Clayton & Hoyle 11974(1 . red giants 
(Norgaard 198fJ, and Wolf- Rayet stars (Dearborn & Blake 
1985 ) . It seems difficult to disentangle the primary galactic 
positron source based only on theoretical grounds, mainly 
due to the (highly) uncertain positron yields, but also due 
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to the uncertain distribution and duty cycle of the source 
populations. 

Help is expected from a detailed study of the 511 keV 
line emission morphology. The celestial 511 keV intensity 
distribution should be tied to the spatial source distri- 
bution, although positron diffusion and effects associated 
with the annihilation physics may to some extent blur this 
link. First estimations of the 511 keV emission morphology 
were obtained by the Oriented Scintillation Spectrometer 
Experiment (OSSE) on-board the Compton Gamma-Ray 
Observatory (CGRO) satellite (Purcell et al. HSMl Cheng 
et al. ITW71 Purcell et al. 119371 Milne et al. I2UUU1 Milne 
et al. I2001f) . but observations were restricted to the inner 
Galaxy, giving only a limited view of the 511 keV emis- 
sion distribution. With the launch of ESA's INTEGRAL 
satellite in October 2002, a new gamma-ray observatory is 
available that allows a detailed study of positron annihi- 
lation signatures. In particular, the imaging spectrometer 
SPI (Vedrenne et al. I2003|) . one of the two prime instru- 
ments on-board INTEGRAL, has been optimised for the 
study of line radiation, combining high-resolution spec- 
troscopy (R ~ 250 at 511 keV) with modest angular res- 
olution (3° FWHM). 

We present in this work an all-sky map of 511 keV 
gamma-ray line emission, with the goals of determining 
the morphology of the emission in the Galaxy and of 
searching for previously unknown sources of 511 keV emis- 
sion anywhere in the sky. The present public data archive 
does not yet cover the entire celestial sphere, but the un- 
exposed regions are limited to a few areas at high galactic 
latitudes, comprising less than 5% of the sky.The resulting 
point-source sensitivity is better than 2x 10 -4 ph cm _2 s _1 
for many regions along the galactic plane, allowing for the 
first time the extraction of information about the distri- 
bution of positron annihilation all over the Galaxy. We 
do not address the distribution of positronium continuum 
emission in this paper, since the subtraction of the dif- 
fuse galactic continuum emission is a distinct data analysis 
challenge. A map of positronium continuum emission will 
be presented elsewhere (Weidenspointner et al., in prepa- 
ration). 

Earlier results on the 511 keV line emission mor- 
phology as observed by SPI have been presented by 
Jean et al. l|2003a[l . Knodlseder et al. (20033), an d 
Weidenspointner et al. I|2004() . and were based on obser- 
vations performed during the galactic centre deep expo- 
sure (GCDE) of 2003. Using a "light bucket" approach 
which neglects the coding properties of the SPI mask, 
Teegarden et al. (2005$ derived upper limits on electron- 
positron annihilation radiation from the galactic disk us- 
ing core-programme data combined with open-programme 
observations at low galactic latitudes (|&| < 20°). In the 
present paper we provide for the first time an all-sky anal- 
ysis using all public data of the first INTEGRAL mission 
year. 

Spectroscopic characteristics of the 511 keV line based 
on SPI data have been published by Jean et al. (2003a), 
Lonjou et al. : 201711) . and Churazov et al. I|20()5[> . We will 



present the 511 keV line profile that we obtain from the 
all-sky dataset elsewhere (Jean et al., in preparation). 

This paper is organised as follows. Section [2] describes 
the observations and the data preparation. Section ex- 
plains the treatment of the instrumental background. In 
section 0] we present the first all-sky map of 511 keV 
gamma-ray line radiation and determine the morphology 
of the emission. Section 0] also describe searches for corre- 
lations with tracers of galactic source populations in order 
to shed light on the origin of the positrons. In sectionOOwe 
discuss the implications of the observations for the galactic 
origin of positrons, and we conclude in section |SJ 

2. Observations and Data Preparation 

The data that were analysed in this work consist of those 
included in the December 10, 2004 public INTEGRAL 
data release (i.e. orbital revolutions 19-76, 79-80, 89-122) 
plus the INTEGRAL Science Working Team data of the 
Vela region observed during revolutions 81-88. The data 
span the IJD epoch 1073.394 - 1383.573, where IJD is the 
Julian Date minus 2 451 544.5 days. 

We screened the data for anomalously high counting 
rates (typically occurring at the beginning and the end of 
an orbital revolution due to the exit and entry of the radia- 
tion belts) and for periods of solar activity (as monitored 
by the SPI anticoincidence system) and excluded these 
periods from the data. This data screening has turned 
out to be crucial for reducing the systematic uncertain- 
ties in the data analysis related to instrumental back- 
ground variations. After data screening, the dataset con- 
sists of 6821 pointed observations, with a total exposure 
time of 15.3 Ms. Typical exposure times per pointing are 
1200 — 3400 seconds, but a few long staring observations 
of up to 113 ks exposure time are also included. 

FigureHshows a map of the resulting effective SPI ex- 
posure at 511 keV. The maximum exposure of 2.1 x 10 8 
cm 2 s occurs towards the galactic centre region thanks 
to data obtained during a long dedicated observation of 
this region. 1 A relatively uniform exposure of ~ 3 x 10 7 
cm 2 s has been achieved for galactic longitudes |Z| < 50° 
and latitudes |6| < 15°. Regions of peculiarly high expo- 
sure (~ 5 x 10 7 cm 2 s) are found in Cygnus, Vela and 
towards the Large Magellanic Cloud. In addition, particu- 
larly well exposed sources (^2x 10 7 cm 2 s) are the Crab 
nebula, 3C 273, NGC 4151, M 94, NGC 936 (during the 
SN2003 gs outburst) and the Coma cluster. Unexposed re- 
gions are found mostly at intermediate galactic latitudes 
(\b\ ~ 30° — 60°), and towards the south galactic pole. 

A map of the resulting narrow-line 3er point-source sen- 
sitivity of SPI at 511 keV is shown in Fig. [3 To evaluate 
the sensitivity, an energy band of 7 keV centred at 511 keV 
has been used. The choice of such a relatively wide band 
eliminates any bias due to the germanium detector degra- 
dation and annealing cycles, as well as any bias/effect due 

1 To obtain the effective exposure time, the exposure has to 
be divided by the effective area at 511 keV of about 75 cm 2 . 
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Fig. 2. SPI narrow line 3<t point-source sensitivity at 
511 keV for the analysis interval 507.5 — 514.5 keV (con- 
tours are labelled in units of 10 -4 ph cm _2 s _1 ). 



to gain calibration uncertainties. It also takes into account 
moderate 511 keV line broadening, as reported by Jean et 
al. I|2003a|l . 

Over large regions of the sky, and in particular in 
the galactic plane, a sensitivity better than 2 x 10 
ph cm s is reached. A best point-source sensitivity of 
5 x 10~ 5 ph cm~ 2 s _1 is achieved towards the galactic cen- 
tre direction. The sensitivity to extended diffuse emission 
becomes slightly worse with increasing emission size, and 
depends on the exposure pattern in the region of interest. 
For example, for a 2d angular Gaussian surface brightness 
distribution centred on the galactic centre, the 511 keV 
line sensitivity worsens from 5 x 10~ 5 ph cm _2 s _1 for a 
galactic centre point-source to 7 x 10~ 5 ph cm _2 s _1 for 
an extended source of 8° FWHM. 



Only single-detector event data have been analysed in 
this work (multiple-detector event data do not contribute 
significantly to the SPI sensitivity at an energy of 511 keV; 
c.f. Roques et al. l2003|) . Energy calibration was performed 
orbit-wise, resulting in a relative (orbit-to-orbit) calibra- 
tion precision of ~ 0.01 keV and an absolute accuracy of 
~ 0.05 keV (Lonjou et al. l2HUl) . 

The data have been analysed by sorting the events in 
a 3-dimensional data-space, spanned by the (calibrated) 
event energy, the detector number, and the SPI pointing 
number. An energy binning of 0.5 keV has been chosen, 
well below the instrumental energy resolution of 2.12 keV 
at 511 keV. 

3. Background modelling 

The most crucial step in SPI data analysis consists of the 
precise modelling of the time variability of the instrumen- 
tal background. In the region of the 511 keV line, the 
instrumental background consists of a nearly flat contin- 
uum and a (broadened) instrumental 511 keV line orig- 
inating from positron annihilation within the telescope 
(Teegarden et al. 2004). Since the time variation of the 
continuum component differs from that of the line com- 
ponent we model them independently. The background 
model for a given data-space bin, indexed in the following 
by the pointing number p, the detector number d and the 
energy bin e, is then given by 

b Pl d,e = b C p°S*e + b p™e i 1 ) 

(note that for the analysis presented in this work a single 
energy bin has been used, covering the energy interval 
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507.5 — 514.5 keV; however, for clarity and reference in 
future works we give here the complete energy-dependent 
formalism) . 

The time variation of the continuum component is ex- 
trapolated from that observed in an continuum energy 
band adjacent to the 511 keV line. We used the energy 
band -E a dj = 523 — 545 keV, situated above the 511 keV 
line, in order to exclude any bias due to positronium con- 
tinuum emission that appears below 511 keV. To reduce 
the statistical uncertainty that arises from the limited 
counting statistics, we smoothed the time variation by lo- 
cally adjusting the rate of saturated events in the germa- 
nium detectors (GEDSAT) to the adjacent counting rate 
(GEDSAT turned out to provide a good first order tracer 
of the background variation in SPI; c.f. Jean et al. 12003b). 
The predicted number of continuum background counts in 
data-space bin (p, d, e) is then given by 

b c °Z = 9 P ,d x T p , d x 

A e l^p'=p-A p 2^e'£E adj n p' ,d,d 

2-e'e£ adJ *V Ep'= P -A p 9 P ',d x T p , A 

where 

— g Pt d is the GEDSAT rate for detector d, averaged over 
the time period spanned by pointing p, given in units 
of counts s _1 , 

— Tp/^d is the lifetime for detector d during pointing p', 
given in units of seconds, 

— A e is the energy bin size for spectral bin e, given in 
units of keV (here A e = 0.5 keV), and 

— n p > ,d, e ' is the number of observed counts for pointing p' , 
detector d, and energy bin e', given in units of counts. 

The number of pointings used for smoothing, given by 
2A p + 1, is determined for each pointing p and detector d 
by satisfying the constraint 

(p+a p \ 
T P'4 ^ T min • (3) 
p'=p-A p J 

An accumulated lifetime of T m j n = 20 hours has shown 
to provide an optimum compromise between reducing 
the statistical uncertainty (due to the limited number of 
events in the adjacent energy band) and reducing the sys- 
tematic uncertainty (due to the fact that the GEDSAT 
rate does not predict the background to infinite precision). 
In other words, continuum background variations shorter 
than ~ 20 hours are modelled by the GEDSAT rate while 
variations on longer time scales are modelled by the ob- 
served event rate in the 523 — 545 keV band. 

The time variation of the line component was modelled 
for each detector d and energy bin e separately using a 
multi-component template of the form 

= C + 0g x 9 P ,d + /3g f 9 d (t')e^'-^dt'. (4) 

J to 



This template consists of a constant term j3\ ' plus the 

GEDSAT rate g p4 scaled by (3 { d 2) e plus the GEDSAT rate 
9d(t') convolved with an exponential decay law, scaled by 
j3 d ' (the convolution integral is taken from the start of the 

INTEGRAL mission t Q up to date t). The coefficients /3 d l) e 
of the template are adjusted during the analysis for each 
SPI detector d and energy bin e using a maximum likeli- 
hood fitting procedure (again, in the analysis presented in 
this paper only a single energy bin is used). The constant 
term has been introduced to provide for non-linearities 
between the background variation and the GEDSAT rate. 
In fact, it turns out that f3 de are negative. An equally 
good background predictor is obtained if the GEDSAT 
rate raised to a power of ~ 1.1 is taken, but using a con- 
stant instead of a powerlaw has the advantage of having 
the background variation template decomposed into a lin- 
ear combination of terms. The third component makes 
provision for a long term build-up that is seen in the in- 
tensity of the 511 keV background line, and that is tenta- 
tively attributed to production of the isotope 65 Zn which 
has a decay time of r = 352 days. The precise value of r is 
in fact weakly constrained by the present data, and a lin- 
ear slope provides an equally good fit of the instrumental 
511 keV line background. 

Although the background model defined by Eqs. QJ- 
J2J, which hereafter is called model DETE, predicts the in- 
strumental background to good accuracy, significant resid- 
uals remain after subtracting off the background model 
and a model of the sky intensity distribution from the 
data (c.f. Fig. We found that these residuals can lead 
to systematic biases in the study of the morphology of 
the 511 keV emission, in particular for the determination 
of the longitude profile of the emission. These biases can 
be explained by the telescope pointing strategy that has 
been adopted for a large fraction of the galactic centre 
deep exposure (GCDE) of the INTEGRAL core program: 
a slow (5°/day) scan of the galactic plane from negative to- 
wards positive longitudes combined with rapid (3° /hour) 
excursions in galactic latitude. As result, the longitude 
profile of the 511 keV line emission is encoded in count 
rate variations on timescales of days while the latitude 
profile is encoded in count rate variations on timescales 
of hours. Within a few hours the SPI instrumental back- 
ground is sufficiently stable to be accurately predicted by 
our model, hence the latitude profile is rather well de- 
termined. However, on timescales of days the background 
variations are more difficult to predict to sufficient accu- 
racy, potentially leading to systematic trends in the deter- 
mination of the longitude profile. 2 

2 In a preliminary analysis in which we treated a much 
smaller dataset, systematic background uncertainties sug- 
gested a significant elongation of the galactic centre bulge emis- 
sion along the galactic plane. This elongation was artificial and 
had been produced by background variations that were not 
fully explained by our model. Removing the short period of 
data with the strongest background variations removed also 
the apparent elongation of the bulge emission. 
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Fig. 3. Residual count rate as function of pointing num- 
ber for the DETE (top) and DRBIT-DETE (bottom) back- 
ground models for the energy band 507.5 — 514.5 keV. 
In addition to the background models the best fitting 2d 
Gaussian surface brightness model (c.f. section 14. 2. 2[) has 
been subtracted from the data. For clarity, the data have 
been rebinned into groups of 50 pointings. The shaded 
area indicates count rate variations of ±0.5%. For com- 
parison, the maximum 511 keV line signal amplitude cor- 
responds to approximately ~ 2% of the observed count 
rate. 



fact that flux measurements require an average determi- 
nation of the count rate level and are not sensitive to the 
count rate contrast. Apparently, the count rate residuals 
approximately average to zero (c.f. Fig. 0). 

We therefore opted for a two step approach where we 
first determine the morphology using ORBIT-DETE, and 
then, using the optimum morphology parameters, deter- 
mine the 511 keV flux using DETE. In this way we re- 
cover the good sensitivity of SPI for 511 keV flux mea- 
surements that was reduced by a factor of ~ 2 by the 
usage of ORBIT-DETE. The comparison of the flux levels 
determined using DETE and ORBIT-DETE provides us with 
a measure of the systematic uncertainty in the flux deter- 
mination, which in general is smaller than the statistical 
uncertainty obtained with DETE. We add the systematic to 
the statistical uncertainty in quadrature and quote the re- 
sult as total error on the flux measurement. In cases where 
uncertainties in the morphology (such as the size of the 
emission region) introduce some uncertainty on the flux, 
we have also added this uncertainty to the total error in 
quadrature. 



In order to improve the background model on long 
timcscales, we studied also a class of models where we ad- 
just the longterm variations during model fitting. For this 

(2) 

purpose we adjust the model parameters f3 d e not only 
for all SPI detectors but also for time intervals of fixed 
duration T. In this way, systematic uncertainties in the 
background model on timescales longer than T are re- 
moved. Fitting the background for each orbital revolution 
(T ~ 3 days) is adequate to reduce systematic trends well 
below the statistical uncertainties (c.f. bottom panel of 
Fig.|nj). This method is similar to the method that we ap- 
plied in our earlier works (Jean et al. EUQSaJ Knodlseder 
et al. l20()4al Weidenspointner et al. l2TY)4l 



with the differ- 
ence that we now also fit the background model for each of 
the SPI detectors separately, and that we included in ad- 
dition a constant term and a build-up term in the model 
(see Eq. 0J. Hereafter this second background model is 
called ORBIT-DETE. 

The introduction of additional parameters in 
ORBIT-DETE with respect to DETE leads to a sub- 
stantial loss in sensitivity. The detection significance of 
galactic centre 511 keV line emission drops from ~ 50er 
for DETE to ~ 22er for ORBIT-DETE. However, it was 
found that the statistical accuracy of the morphology 
determination, which is driven by the count rate contrast 
in the data-space rather than the count rate level, is not 
degraded by the introduction of additional parameters, as 
long as T > 2 days. Consequently, using the ORBIT-DETE 
model for the morphological characterisation of the 
511 keV line emission is the optimum choice that keeps 
a high statistical accuracy while reducing the systematic 
uncertainties in the analysis. 

On the other hand, despite the systematic uncertain- 
ties, DETE is accurate enough to allow for a precise deter- 
mination of 511 keV line flux levels. This is related to the 



4. Results 

4.1. Imaging 

To determine a model independent map of the 511 keV 
gamma-ray line intensity distribution over the sky, we em- 
ployed the Richardson-Lucy algorithm (Richardson 119721 
Lucv llUT^jl . This type of algorithm is widely used for im- 
age deconvolution, and has in particular been successfully 
employed for the analysis of gamma-ray data of CGRO 
(Knodlseder et al. U999al Milne et al. 12710711 . 

We implemented the accelerated version ML-LINB-1 
of Kaufman (1987) of the Richardson-Lucy algorithm for 
our analysis, which iteratively updates the sky intensity 
distribution /* 



/• +1 using the relation 



f- +1 = f- + '\' ir jfj 



'x^ N ( m. 



1 )Ri 



(5) 



where Rij is the instrumental response matrix (linking the 
data space, indexed by i, to the image space, indexed by 
j), rii is the number of counts measured in data space 
bin i, e\ = J2j=i R ijfj + i s tne predicted number of 
counts in data space bin i after iteration k (pi being the 
predicted number of instrumental background counts for 
bin i), N and M are the dimensions of the data and image 
space, respectively, and A fc is an acceleration factor that is 
obtained by constrained maximum likelihood fitting (with 
the constraint that the resulting sky intensities remain 
positive) . 

To avoid noise artefacts in the weakly exposed regions 
of the sky, we weighted the image increment with a quan- 
tity that is related to the sensitivity of the instrument, 
given by Wj — (Yli-i Rij) 1 ^ 2 - We verified that introduc- 
ing this weighting had no impact on the image reconstruc- 
tion in the well exposed regions of the sky. In addition, we 
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Fig. 4. Richardson-Lucy image of 511 keV gamma-ray line emission (iteration 17). Contour levels indicate intensity 
levels of 10~ 2 , 10 -3 , and 10~ 4 ph cm _2 s _1 sr _1 (from the centre outwards). 
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Galactic longitude (deg) Galactic latitude (deg) 

Fig. 5. Longitude and latitude profiles of the image shown in Fig. 0] (integration range \l\ < 30°, |6| < 30°). 



smoothed the iterative corrections on the right hand side 
of Eq. fusing a 5° x 5° boxcar average. In this way the ef- 
fective number of free parameters in the reconstruction is 
reduced and image noise is damped to an acceptable level. 
The application of more sophisticated image reconstruc- 
tion methods involving wavelet based multi-resolution al- 
gorithms aiming at a complete suppression of image noise 
(Knodlseder et al. 1999a) will be presented elsewhere. 

The resulting all-sky image of the 511 keV line emis- 
sion is shown in Fig. ^ longitude and latitude profiles 
of the emission are shown in Fig. We have chosen to 
stop the iterative procedure after iteration 17 since at this 
point the recovered flux and the fit quality correspond 
approximately to the values that we achieve by fitting as- 
trophysical models to the data (c.f. section FOll . In this 
way we make sure that we are not in the regime of over- 
fitting, which is characterised by substantial image noise 



and artificial image structures. On the other hand, sim- 
ulations showed that faint diffuse emission, as expected 
for example for a galactic disk component, would not be 
recovered at this point. 

Figure 01 reveals that the 511 keV sky is dominated by 
prominent emission from the bulge region of the Galaxy. 
Beyond the galactic bulge, no additional 511 keV emis- 
sion is seen all over the sky, despite the good exposure in 
some regions (e.g. Cygnus, Vela, LMC, anticentre, north 
galactic pole region). The 511 keV emission appears sym- 
metric and centred on the galactic centre, with indica- 
tions for a slight latitude flattening. The latitude flatten- 
ing could be either due to an inherent asymmetry of the 
bulge component or due to the presence of an underlying 
faint galactic disk component. Indeed, if the Richardson- 
Lucy iterations are continued, a faint disk-like structure 
emerges (c.f. Fig. - Yet the image starts to become pol- 
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luted by noise and we cannot exclude the possibility that 
the apparent disk emission is artificially created by the ex- 
posure pattern that follows the galactic plane. Therefore 
we employ more quantitative methods in the next section 
to assess the significance of the possible disk emission. 

By fitting Gaussian functions to the longitude and lat- 
itude profiles of the image (c.f. Fig. |3J) we estimate the 
extent of the emission to ~ 13° x 10° (FWHM). Figure 
indicates, however, that the emission profiles are not well 
represented by Gaussian functions. The emission is better 
described by a compact (FWHM ^5°) core and a more 
extended halo (FWHM 10° - 20°). We want to empha- 
sise, however, that this qualitative analysis should not be 
pushed too far, since image deconvolution is a non-linear 
process which is easily affected by image noise and expo- 
sure biases. 

4.2. Morphologies! characterisation 
4.2.1. Method 

To make a quantitative assessment of the morphology 
of the 511 keV line emission we use a maximum likeli- 
hood multi-component model fitting algorithm. Assuming 
Poisson noise for the measured number rii of events in each 
of the N data-space bins, the algorithm maximises the log 
likelihood 

JV 

In L = V^Jtj In e.% — — In nil, (6) 

i=l 

where = J2k akS i + ^iiP) ^ s the predicted number 
of (source plus background) counts in data space bin i, 
s!- = J2j=i fjRij i s tnc sky intensity model /* folded 
into the data space (Rij being the instrumental response 
matrix), bi(f3) is the background model (c.f. Fig. HJ, and 
ah and (3 are scaling factors for the sky intensity and the 
background model, respectively, that are adjusted by the 
fit. 

Detection significances (and parameter errors) are 
estimated using the maximum likelihood ratio test 
fCash rH)79[l . We calculate the maximum log likelihood- 
ratio MLR = — 2(lnL —hxLi) between two models (hy- 
potheses) , where for the first one we constrain a number p 
of the parameters to specific values (resulting in Lq) while 
for the second one all parameters are left free (resulting in 
Li). In the case that L\ provides a satisfactory fit of the 
data, MLR is then distributed like a \p distribution with p 
degrees of freedom. Statistical parameter errors were esti- 
mated using the formalism of Strong (1985). Throughout 
this paper the error bars quoted are la. 

We call the maximum log likelihood-ratio (MLR) of a 
model the difference between the log likelihood obtained 
by fitting all model parameters and the log likelihood ob- 
tained by fitting only the background model to the data 
(i.e. for L\ all parameters ak and (3 vary freely while for 
Lq all ak are constrained to zero and only the f3 are al- 
lowed to vary). To compare models with different num- 
bers of free parameters, we quote the reduced maximum 



Table 1. Morphology of the emission assuming a 2d an- 
gular Gaussian surface brightness distribution. 



Quantity Measured Value 



RMLR (DOF) 462.2 (5) 

lo -0.6°±0.3 C 

b +0.1°±0.3 C 

Al (FWHM) 8.1° ±0.9° 

Ab (FWHM) 7.2° ± 0.9° 

511 keV flux (1(T 3 ph cm'V 1 ) 1.09 ± 0.04 



log likelihood-ratio, RMLR = MLR- DOF, with DOF be- 
ing the number of free parameters ak of the sky intensity 
model. 

4.2.2. 2d surface brightness distribution 

As a first step we characterise the apparent morphology 
of the 511 keV line emission on the sky using a 2d angu- 
lar Gaussian surface brightness distribution for which we 
determined the centroid, lo,bo, the longitude and latitude 
extent, AZ,A6, and the 511 keV line flux. The results of 
this analysis are summarised in Table ^ the best fitting 
model intensity distribution is shown in Fig. |H1 

The analysis confirms our earlier findings (Jean et 
al. I2003al Knodlseder et al. I20()4al Weidenspointner et 
al. I2()()4|) of a compact and symmetric 511 keV line emis- 
sion distribution towards the galactic centre. The centroid 
of the emission appears slightly offset from the galactic 
centre direction, at the statistical 2er level, but we do 
not claim that this offset is significant. From our earlier 
analyses we learned that the centroid can be shifted by 
this amount simply from the combined effect of statistical 
and systematic biases in the modelling of the instrumental 
background. 

Within the statistical uncertainties, the emission ap- 
pears fully symmetric, with an extension of ~ 8° 
(FWHM). Formally, we determine a marginal emission 
flattening of Ab/Al = 0.89 ± 0.14. The total 511 keV flux 
is (1.09±0.04) x 10~ 3 ph cm _2 s -1 , where the quoted error 
includes the uncertainty in the extent of the emission and 
the statistical and systematic measurement errors (c.f. sec- 
tion EJ). 

The RMLR of 462.2 that has been obtained using the 
0RBIT-DETE background model converts into a formal de- 
tection significance of 22c. Using the DETE background 
model and including the systematic uncertainties results 
in a substantially higher detection significance of 34cr. 
Neglecting systematic uncertainties would even boost the 
detection significance towards 49er. 

4.2.3. Galactic models 

To determine the galactic positron-electron annihilation 
rate requires modelling the spatial distribution of the 
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positron-electron annihilation. The 511 keV photon lu- 
minosity L511 is related to the positron luminosity L p 
through L511 = (2 — 1.5/ p ) x L p where / p is the positron- 
ium (Ps) fraction, defined as the fraction of positrons that 
decay via positronium formation (Brown, & Leventhal 
HW7\ . Using / p = 0.93 ± 0.04 that has been determined 
from OSSE observations (Kinzer et al. 12001(1 results in a 
conversion from 511 keV photon luminosity to a positron- 
electron annihilation rate of L p = (1.64 ± 0.06) x L511. 

We here compare models of bulge, disk, and halo com- 
ponents with the data. Based on galactic model density 
distributions p(x, y, z) we calculate the expected all-sky 
511 keV intensity f(l,b) towards direction (l,b) by inte- 
grating the volume emissivity p(x, y, z) along the line of 
sight s: 

f(l,b) = ^J p(x,y,z)ds (7) 

(the galactic centre has been assumed to be at a distance 
of Rq — 8.5 kpc). Galactic 511 keV photon luminosities 
are calculated by integrating p(x,y,z) over the galactic 
volume, 

L p = J p(x,y,z)s 2 dsdn (8) 

assuming an outer Galaxy radius of i? m ax = 15 kpc. 

Since the 511 keV line emission is primarily arising 
from the galactic centre region we fitted in a first step 
models of the galactic stellar bulge to the data. To account 
for uncertainties in our knowledge about the morphology 
of this component (which are related to our location in the 
galactic plane amid the obscuration by interstellar dust) 
we compared a variety of proposed bulge models to the 



data. The models were gathered from Dwck et al. (1995) 
and Freudenreich ( 1998), who modelled the distribution of 
K and M giant stars using DIRBE near-infrared skymaps, 
and from Picaud & Robin fZQMfl who analysed data from 
the DENIS near-infrared survey. There is an accumulat- 
ing body of evidence that the stellar distribution in the 
bulge is bar-shaped, and except for models GO and E0, all 
employed bulge models have triaxial morphologies that 
differ in the orientation angles, the scale lengths, and the 
radial density profiles. Details of the models are given in 
Appendix^ the results of the analysis are summarised in 
Table |2 and best fitting 511 keV intensity distributions 
are shown in Fig. |S| 

The best fitting bulge models are E3, G3, Sf, and Ep. 
Reasonably good fits are also obtained for Pp and Spr, 
while only moderate fits are achieved for the remaining 
models. Our ranking is similar to that established from 
the analysis of the DIRBE and DENIS near-infrared data 
(Dwek et al. Freudenreich ITM81 Picaud & Robin 

2004). The best fitting bulge models fit the data as well 
as the adjusted 2d angular Gaussian surface brightness 
distribution. This means that models of the galactic stellar 
bulge are able to explain satisfactorily the morphology of 
the 511 keV bulge emission. 

In a second step we fitted the 511 keV emission using 
parametric models of the galactic bulge and halo mor- 
phology in order to determine the scale of the emission. 
For the bulge models GO' and E0' we adjust the radial 
scale length (Rq) and vertical scale height (zq) 1 while for 
the galactic halo model H' we determine the density slope 
powerlaw index (n), the inner cutoff radius (a c ), and the 
axis ratio (e). In addition, we employed a model composed 
of a set of galactocentric nested shells of constant density 
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Table 2. Galaxy model fitting results (see text). The columns give (1) the model, (2) the RMLR (obtained using 
the ORBIT-DETE background model) and the number of free model parameters (DOF), (3) the bulge scale length, 
(4) the bulge scale height, (5)-(8) the 511 keV photon luminosity of the model components, and (9)-(12) the total 
47T integrated 511 keV line all-sky flux in each of the model components. The prime indicates model components for 
which the scaling parameters were adjusted by the fit. The figures in parenthesis quoted in columns (5)-(12) indicate 
la uncertainties in the last digit. 
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(model 'Shells') to determine the radial density profile of 
the 511 keV emission. We varied the radii of the shells and 
the number of shells in order to maximise the MLR, whilst 
limiting the number of shells to the minimum required to 
satisfactorily describe the data. 

The data suggest a symmetric bulge emission profile, 
with scale lengths between 300 and 600 pc. The RMLRs 
are comparable to the best fitting bulge models that we 
tested before. The data are equally well fitted by a model 
of the galactic halo, with a density powerlaw index of n = 
3.0 ± 0.3, an inner cutoff radius of a c — 0.39 ± 0.08 kpc, 
and a flattening of e = 0.81 ± 0.12. Most studies of the 
stellar halo population suggest power indices between 2.4 
and 3.5 and flattenings in the range 0.6 to 1.0, while the 
inner cutoff radius is basically undetermined (Robin et 
al-ED00 and references therein). Our values are compatible 
with those of the stellar halo population, but the large 
uncertainties in the stellar halo morphology do not allow 
firm conclusions to be drawn. 

The nested shell model provides the best fit to the data 
thanks to its flexibility in adjusting the radial density pro- 



file of the emission. A satisfactory fit is achieved by using 
two shells with radii — 0.5 and 0.5 — 1.5 kpc; splitting up 
these shells in a finer binning, moving the shell interface 
radius or adding more shells does not significantly improve 
the fit. In particular, we detect no significant 511 keV 
bulge emission from galactocentric distances 1.5 kpc. 
The radial dependence of the 511 keV volume emissivity 
is plotted in Fig. [7] For illustration we added the result 
of a third shell to the figure that covers radial distances 
of 1.5 — 3.0 kpc and for which the flux is consistent with 
zero. Our fit reveals a drop in the annihilation emissivity 
by one order of magnitude between the inner — 0.5 kpc 
and the outer 0.5 — 1.5 kpc shell, confirming the existence 
of a narrow core plus an extended halo of 511 keV emission 
that has already been suggested by the imaging analysis 
(c.f. section EtJ). 

In a third step we added galactic disk components to 
the bulge and halo models. For the galactic disk we tested 
models of young (model DO) and old (model Dl) stellar 
populations (Robin et al. 2003 ). With both models we find 
clear evidence for 511 keV line emission from the galac- 
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Fig. 7. Radial dependence of the 511 keV volume emis- 
sivity as derived from the bulge model 'Shells'. 

Table 3. Summary of model fitting results. Fluxes are 
given as total 4tt integrated all-sky values. Annihilation 
rates have been calculated assuming / p = 0.93. 



Quantity 


Bulge 


Halo 


Disk 


Flux (10~ 3 ph cni'V 1 


1.05 ±0.06 


1.6 ±0.5 


0.7 ±0.4 


L511 (10 43 ph s" 1 ) 


0.90 ± 0.06 


1.2 ±0.3 


0.2 ±0.1 


L p (10 43 s- 1 ) 


1.50 ±0.10 


2.0 ±0.5 


0.3 ±0.2 


Total flux (10~ 3 ph cm- 
Total L 5 n (10 43 ph s _1 ; 
Total L p (10 43 s" 1 ) 
B/D flux ratio 
B/D luminosity ratio 


V 1 ) 




1.5- 2.9 
1.0 - 1.7 

1.6- 2.8 
1-3 

3-9 



most equally good fit to the data. Due to their degeneracy 
fitting both simultaneously is not meaningful. 

The halo model leads to a considerably larger flux, lu- 
minosity and annihilation rate than the bulge model due 
to the presence of a flat and extended tail in this distribu- 
tion (c.f. Fig.[SJ. Currently, our data do not allow to detect 
this tail, and thus, they do not allow to discriminate be- 
tween bulge and halo models. Future deep observations 
at intermediate galactic latitudes that are scheduled for 
the INTEGRAL AO-3 observing period aim in measuring 
this emission tail, promising to provide constraints that 
will allow in the future to disentangle between the differ- 
ent emission morphologies. 

The data suggest bulge-to-disk 511 keV flux ratios in 
the range 1 — 3, where the lower boundary is obtained 
for the short scale-length old stellar disk model Dl which 
suggests larger disk flux values than the young stellar disk 
model DO. Halo-to-disk 511 keV flux ratios are even larger, 
in the range 2 — 4, owing to the larger flux in the halo com- 
ponent. The large uncertainty in these ratios arises from 
the low intensity of the galactic disk component, which for 
the analysed dataset is just above the SPI detection limit. 

We also note that the bulge-to-disk 511 keV photon 
luminosity ratio is much higher than the bulge-to-disk 
flux ratio and lies in the range 3 — 9. This difference is 
explained by the fact that the average squared distance 
s 2 = f p(s)s 2 dsdS}/ J p(s)dsdtt, which defines the dis- 
tance at which a source of luminosity L p produces the 
observed 511 keV line flux, is smaller for the galactic disk 
than for the galactic bulge. In other words, to produce the 
same 511 keV flux at Earth, the intrinsic luminosity of the 
bulge has to be larger than that of the disk. 4 It is there- 
fore important to quote explicitly the quantity for which 
we discuss the bulge-to-disk ratio. The same rational also 
holds for the halo-to-disk 511 keV photon luminosity ra- 
tio, which is larger than the corresponding flux ratio, and 
which is comprised between 4 — 13. 



tic disk. Adding disk models DO and Dl to bulge or halo 
models consistently improves the fit leading to a detection 
of the disk emission at the 3 — 4cr level. 3 Formally, Dl pro- 
vides a better fit than DO, but the difference is marginal. 
The signal from the disk is still too faint in our present 
dataset to deduce anything about its morphology. 

The flux, luminosity and annihilation rate in the bulge, 
halo and disk components are summarised in Table |31 
Recall that the bulge and halo components are alterna- 
tives and their contribution should not be added to derive 
the total galactic values. Either component provides an al- 

3 Using the quoted RMLRs, the formal significance of the 
disk emission amounts only to 2 — 3cr. However, the tabulated 
RMLRs have been obtained using the ORBIT-DETE background 
model which is less sensitive to 511 keV line emission than 
the DETE model. Using the procedure outlined in section [3] we 
reduce the flux uncertainties and increase the detection signif- 
icance to 3 — 4a. 



4.3. Correlation with tracer maps 

To gain insight into the nature of the galactic positron 
sources, we searched for correlations between the 511 keV 
line emission morphology and all-sky intensity distribu- 
tions observed at other wavelengths. This work was in- 
spired by a similar study that Knodlseder et al. (1999b) 
performed to understand the morphology of the 1.8 MeV 
gamma-ray line emission (arising from radioactive decay 
of 26 A1) observed by the COMPTEL telescope aboard 
CGRO. Through their analysis, the authors could estab- 
lish a tight correlation between the morphology of galac- 
tic microwave free- free emission and that of 1.8 MeV line 
emission, hinting towards a massive star origin of 26 A1. 

4 The difference between bulge-to-disk flux and luminosity 
ratio is only important for our home Galaxy and is related to 
the fact that the Sun is located within the galactic radius. For 
external galaxies this difference disappears since their bulge 
and disk appear at the same distance to us. 







Fig. 8. All-sky maps of the best fitting models of 511 keV gamma-ray line emission (see text for a description of the 
models). Contour levels indicate intensity levels (from the centre outwards) of 10 -2 , 1CP 3 , and 1CP 4 ph cm~ 2 s _:1, sr . 
The resulting RMLRs of the model fits are quoted in the upper-right corner of the panels. 



The tracer maps used for the comparison are those listed 
in Knodlseder et al. (1999b). For a detailed description of 
the maps the reader is referred to that work. 

FigureElsummarises the result of the correlation study 
for our 511 keV dataset, where we show the RMLR as a 
function of the tracer map (ordered by increasing photon 
energy or decreasing wavelength) . None of the tracer maps 
is consistent with the data. The maximum RMLR that 
is reached (353.7 for the DIRBE 2.2 /im map) is more 
than 100 units smaller than the values obtained for the 
parametric models of the previous sections. Apparently, 



the 511 keV emission morphology is unique and cannot be 
represented by any known celestial intensity distribution. 

Nevertheless, Fig.[5]shows a clear trend, where the data 
favour maps in the near-infrared domain (DIRBE 1.25 fim 
- 4.9 /im) and the hard X-ray band (HEAO-1) over maps 
observed at longer wavelengths. In particular, the worst 
fits are obtained in the microwave and far-infrared do- 
main where the skymaps trace the young stellar popula- 
tion, either through their ionising radiation (DMR maps 
at v 53 GHz), or through their related molecular gas 
(CO) and cold dust emission (DIRBE 100 \im - 240 /jm). 
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Fig. 9. Reduced maximum log likelihood-ratio (RMLR) 
as function of the tracer map, ordered by increasing pho- 
ton energy (or decreasing wavelength). The dotted line 
illustrates the impact of zodiacal light contamination in 
the DIRBE 12 /im and 25 /im skymaps, the dashed line 
indicates maps that are heavily affected by galactic ab- 
sorption. 



From this it is clear that the bulk of the 511 keV emission 
is not related to the young massive stellar population of 
the Galaxy. 

On the contrary, all best fitting tracers maps show 
the characteristic features of an old stellar population: a 
strong bulge component combined with a short scale ra- 
dius disk component. Apparently, the 511 keV data tend 
to favour such morphologies. This is illustrated in Fig. I1UI 
where we plot the RMLR as function of the bulge-to-disk 
flux ratio of the tracer map, defined as the flux contained 
within a circular region of 6° in radius around the galac- 
tic centre, divided by the flux within galactic latitudes 
b = ±20° outside the circular bulge region (note that the 
precise B /D value depends of course on the exact defini- 
tions chosen for the two regions, but our purpose here is 
to illustrate a trend). Clearly, there is a strong correla- 
tion between the B/D flux ratio and the RMLR, in the 
sense that the larger the B/D flux ratio, the larger the 
RMLR. In particular, the DIRBE near-infrared maps and 
the HEAO-1 map show the largest B/D flux ratios of all 
tracer maps (~ 0.2). Thus, finding a population of objects 
which show a large B/D ratio could provide the key for 
finding the galactic source of positrons. 

4.4. Point-source search 

The modest angular resolution of SPI of about 3° 
(Vedrenne et al. [5003 ) makes it difficult to distinguish 
between point source, point-like, and small-scale diffuse 
emission. So in principle we cannot exclude the possibility 
that the 511 keV gamma-ray line emission that is seen 
towards the galactic bulge region is made of a limited 
number of point sources that blend to simulate diffuse 
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Fig. 10. Reduced maximum log likelihood-ratio (RMLR) 
as function of the bulge-to-disk flux ratio of the tracer 
map. The results for the parametric bulge+disk models 
are also shown. 



emission. In due course the question of any point source 
contribution to the flux will be best addressed using data 
from the imager IBIS on INTEGRAL, which is relatively 
insensitive to diffuse emission, in conjunction with that 
from SPI. Such work is underway and will be reported 
on separately. We here limit ourselves to the constraints 
which can be placed from SPI data alone on such a con- 
tribution. 

We have therefore used the SPIROS algorithm 
(Skinner & Connell 2003) to search our dataset for the 
positions and fluxes of point sources that are compatible 
with our data. SPIROS searches for the most probable po- 
sition for a point source and fits a source at that position 
before repeating iteratively the search using the residu- 
als after sources already found are taken into account. At 
each iteration the positions and fluxes of all sources that 
have been found are optimised by maximising a goodness 
of fit parameter (the x 2 statistic was used here). 

If the 511 keV emission is intrinsically diffuse then ap- 
plication of this algorithm will lead to sources being placed 
at selected positions (regions of high flux and local noise 
peaks) until a distribution of emission is found that is 
consistent (given limited statistics) with the data. Such a 
description is unlikely, however, to be unique and if most 
or all of the flux is from diffuse emission then the par- 
ticular source positions found will have no astrophysical 
significance. We therefore do not present here the detailed 
results of this blind source search and we restrict our- 
selves to discussion of the general conclusions which can 
be drawn from the analysis. 

For models with 7 point sources based on the 
SPIROS solution after iteration 7 we obtain RMLR = 
462.5, slightly inferior to our best fitting diffuse models. 
Subsequent iterations suggest point sources at the edge of 
the exposed regions which are obviously spurious. Only af- 
ter iteration 13 is another point source found in the galac- 
tic bulge region. Fitting this source together with the 7 
sources found earlier leads to RMLR = 471.9, comparable 
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to the best fitting diffuse models. The total flux attributed 
to the 8 point sources is 1.1 x 10~ 3 ph cm -2 s -1 , compara- 
ble with the values obtained for the bulge component using 
our best fitting diffuse models (c.f. Table We therefore 
simply conclude that at least 8 point sources would be 
needed to satisfactorily describe the SPI data. 5 

In addition to the blind search for point sources, we 
also looked for evidence of 511 keV gamma-ray line emis- 
sion from a list of potential candidate objects. Our list 
comprises compact objects, pulsars, supernova remnants, 
star forming regions, globular clusters, nearby (active) 
galaxies, and galaxy clusters. Depending on the expected 
source extent, we searched for either point source emission 
or extended emission, modelled by a 2d angular Gaussian 
surface brightness distribution for which we specify the 
centroid and the FWHM extension. 

The results of the analysis are summarised in Table 
0J None of the sources we searched for showed a signifi- 
cant 511 keV flux, hence we only quote (3<r) upper limits 
in Table 01 Since the emission of the Galaxy may inter- 
fere with the emission from the specific sources (due to 
the large field of view of the SPI instrument), we also in- 
cluded models for the diffuse galactic 511 keV emission for 
the source search. We have used combinations Shells+DO, 
Shclls+Dl, H'+DO and H'+Dl to cover the range of plau- 
sible best fitting diffuse models (c.f. Fig. |SJ and quote al- 
ways the most conservative flux limit. The (less sensitive) 
DRBIT-DETE background model has been used to ensure 
that systematic uncertainties are negligible. 

5. Discussion 

5.1. Comparison with earlier measurements 

Observations of the 511 keV line emission have been made 
by a large number of balloon and satellite borne tele- 
scopes, yet only a few of them provided constraining in- 
formation on the emission morphology. The OSSE instru- 
ment that flew during 1991-2000 on-board CGRO accu- 
mulated so far the largest database for studying the 511 
keV line intensity distribution. The observations of the 
Gamma-Ray Spectrometer on-board the Solar Maximum 
Mission (SMM) (1980-1988) and of the Transient Gamma- 
Ray Spectrometer (TGRS) on-board the WIND mission 
(1995-1997) have also been used to estimate the overall 
511 keV line flux and maximum emission size (Purcell et 
al. Kinzer et al. Tueller et al. Cheng et 

al. ITMTl Purcell et al. ITMTl Harris et al. IT53SI Milne et 
al. l2UUUl Kinzer et al. l2UUT|> . 

The picture that emerged prior to the INTEGRAL 
launch was the following. From the general trend that in- 
struments with larger fields-of-view show larger fluxes it 
was inferred that the 511 keV emission is extended. OSSE 
observations strongly exclude a single point-source located 
at the GC (Purcell et al. lT5Mj> . The OSSE observations 

5 Given the angular resolution of SPI, compact sources with 
an extent ~ 1° — 2° would be considered as point sources in 
this context. 



Table 5. Comparison of SPI results with OSSE measure- 
ments. The bulge parameters £o> bo, Al, and Ab of OSSE 
were taken from Kinzer et al. H2001fl . 



Quantity 


SPI 


OSSE 


lo 


-0.6° ±0.3° 


-0.25° ± 0.25° 


bo 


+0.1° ± 0.3° 


-0.3° ± 0.2° 


Al (FWHM) 


8.1° ± 0.9° 


6.3° ± 1.5° 


Ab (FWHM) 


7.2° ± 0.9° 


4.9° ± 0.7° 


Flux (1(T 3 ph cm" V 1 ) 


1.5 - 2.9 


1 - 3 


B/D flux ratio 


1 - 3 


0.2 - 3.3 



suggest at least two emission components, one being a 
spheroidal bulge and the other being a galactic disk com- 
ponent. A third component, named the Positive Latitude 
Enhancement (PLE), situated about 9° — 12° above the 
GC has been reported (Cheng et al. 119971 Purcell et 
al. ^997), but the morphology and intensity of this com- 
ponent was in fact only poorly determined by the data 
(Milne et al. 12009 • 

The emission, which showed no significant offset from 
the GC, was well fitted by either a model comprising 
a narrow (5 — 6° FWHM) Gaussian bulge plus ~ 35° 
FWHM Gaussian and CO-like disk components, or by 
a centre-truncated R 1 / 4 spheroid plus exponential disk 
model (Purcell et al. HWH Milne et al. l27)(Y)l Kinzer et 
al. l2001[) . The total 511 keV gamma-ray line flux was esti- 
mated to be (1 - 3) x 10~ 3 ph cm s . The distribution 
of flux between the bulge and disk components was only 
weakly constrained by the observations, and depended 
sensitively on the assumed bulge shape. In particular, es- 
timates for the bulge-to-disk (B/D) flux ratio varied from 
0.2 — 3.3 depending upon whether the bulge component 
features a halo (which leads to a large B/D ratio) or not 
(Milne et al. l2lTO|l . 

Our analysis basically confirms the pre-INTEGRAL 
observations (c.f. Table|SJ). One difference is that the bulge 
appears slightly larger in our analysis when compared to 
the OSSE result. We note that OSSE performed differen- 
tial measurements using its 4° x 11° collimator which may 
bias the results towards small values (Kinzer et al.[2001), 
but in any case, the discrepancy is not very significant and 
is not surprising in view of possible systematic uncertain- 
ties. 

Another difference with respect to OSSE is that we 
find no evidence for a feature resembling the PLE. Fitting 
a model of the PLE (2d Gaussian of 5° FWHM located 
at I = —2° and b = 8°) on top of the bulge results in 
a 3ct upper flux limit of 1.5 x 10~ 4 ph cm~ 2 s _1 for the 
PLE. The OSSE team has gradually reduced their esti- 
mates of the flux and significance attributed to this emis- 
sion feature from 5x 10 -4 ph cm~ 2 s _1 (Purcell et al. H997[) 
down to an upper limit of 1 x 10 -4 ph cm _2 s _1 (Milne et 
al. l2001J) . Recently it has been suggested that data analy- 
sis problems linked with variable continuum emission may 
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Table 4. 511 keV narrow line 3tr upper flux limits for selected potential positron sources. 



Source name 


1 


b 


FWHM 


3<7 flux limit 


Type of object 




(dee) 






(10~ 4 ph cm" 2 s _1 ) 




Sgr A* 


0.00° 


+0.00° 




< 1.0 


BHC 


GRS 1758-258 


4.51° 


-1.36° 


_ 


< 0.7 


AtQSO 


Cyg X-l 


71.33° 


+3.07° 


_ 


< 1.0 


BHC 


LS I+61°303 


135.68° 


+1.09° 


_ 


< 3.3 


AtQSO 


IE 1740.7-2942 


359.15° 


-0.12° 


_ 


< 0.9 


MQSO 


Vela X-l 


263.06° 


+3.93° 


- 


< 1.1 


HMXB 


GX 5-1 


5.08° 


-1.02° 




< 0.7 


LMXB 


GRS 1915+105 


45.37° 


-0.22° 


- 


< 1.0 


LMXB 


A 0620-00 


209.96° 


-6.54° 


- 


< 3.8 


LMXB 


Nova Muscae 


295.30° 


-7.07° 


- 


< 2.0 


LMXB 


Cir X-l 


322.12° 


+0.04° 




< 1.1 


LMXB 


Cen X-4 


332.24° 


+23.89° 




< 1.7 


LMXB 


GX 349+2 


349.10° 


+2.75° 




< 0.8 


LMXB 


Sco X-l 


359.09° 


+23.78° 


_ 


< 1.5 


LMXB 


Crab 


184.56° 


-5.78° 


- 


< 1.3 


Pulsar 


Geminga 


195.13° 


+4.27° 


- 


< 2.2 


Pulsar 


PSR J0737-3039 


245.24° 


-4.50° 




< 3.0 


Pulsar 


Kepler 


4.52° 


+6.82° 


- 


< 0.7 


SNR 


Cas A 


111.73° 


-2.13° 


- 


< 1.9 


SNR 


Tycho 


120.07° 


+ 1.44° 




< 2.0 


SNR 


SN 1987A 


279.70° 


-31.94° 




< 1.0 


SNR 


SN 1006 


327.58° 


+14.59° 




< 1.3 


SNR 


Lupus Loop 


329.80° 


+16.00° 


- 


< 1.3 


SNR 


Cygnus 


79° 


0° 


5° 


< 1.8 


Star forming region ( Al source) 


Vela 


265° 


0° 


5° 


< 2.6 


Star forming region ( Al source) 


Carina 


286.5° 


+0.5° 




< 2.2 


Star forming region ( 2o Al source) 


M 22 


9.89° 


-7.55° 


- 


< 0.8 


Globular Cluster 


Palomar 13 


87.10° 


-42.70° 




< 7.0 


Globular Cluster 


uj Cen 


309.10° 


+14.97° 




< 1.7 


Globular Cluster 


NGC 6397 


338.16° 


-11.96° 




< 1.2 


Globular Cluster 


M 4 


350.97° 


+15.97° 


_ 


< 1.1 


Globular Cluster 


Sgr dwarf 


5.61° 


-14.10° 


8° 


< 1.7 


Dwarf Galaxy 


M31 


121.17° 


-21.57° 


3° 


< 11.8 


Galaxy 


LMC 


280.47° 


-32.89° 


10.8° 


< 4.7 


Irregular Galaxy 


QSO B2251-179 


46.20° 


-61.33° 




< 2.9 


QSO 


3C 273 


289.95° 


+64.36° 




< 1.2 


QSO 


3C 279 


305.10° 


+57.06° 




< 1.3 


Blazar 


Cen A 


309.52° 


+ 19.42° 




< 1.7 


Radio Galaxy 


Coma cluster 


58.08° 


+87.96° 


4° 


< 2.0 


Galaxy Cluster 


Perseus cluster 


150.58° 


-13.26° 


4° 


< 4.5 


Galaxy Cluster 


Virgo cluster 


281.63° 


+75.18° 


12° 


< 2.5 


Galaxy Cluster 



account for the reported PLE (Milne 2004), so perhaps our 
non-detection of a PLE feature is not surprising. 

Until now, there have been very few published upper 
limits on 511 keV gamma-ray line emission from point 
sources which take account of diffuse emission. Examples 
are the 3a limits of 1.6 x 10~ 4 ph cm^s" 1 for IE 1740.7- 
2942 (Purcell et al. lTTM| and - 1.4 x 10~ 4 ph cm^s -1 for 



inner Galaxy sources (Milne et al. 2001). Our upper limits, 
summarised in Table 0] are somewhat more stringent. 

Finally, we want to mention that the method of anal- 
ysis used in this work assumes that the 511 keV line 
emission is not time variable. From our analysis of SPI 
data alone we have no indications for time variability. 
Furthermore, OSSE and TGRS measurements revealed no 
significant time variability (Purcell et al. 119971 Harris et 
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al. I1998|l and in addition, our 511 keV line flux measure- 
ments are consistent with those of OSSE and TGRS. Thus 
we believe that our assumption of non- variable 511 keV 
line emission is reasonable. We have, however, not yet per- 
formed a thorough analysis on all relevant timescales. 

5.2. General considerations 

The most distinctive morphological feature of the 511 keV 
emission is the large B/D luminosity ratio of 3 — 9. Unless 
there is a mechanism that strongly suppresses positron 
annihilation in the galactic disk, or that somehow trans- 
ports positrons from the disk into the galactic bulge or 
halo where they annihilate, the positron source popula- 
tion we are seeking for should also exhibit such a high 
B/D ratio. 

The B/D luminosity ratio of 3 — 9 is considerably 
larger than the B/D mass ratio of 0.3 — 1.0 of our 
Galaxy (e.g. Caldwell & Ostriker HMTl Freudenreich l 19981 
Bissantz & Gerhard 12177121 Robin et al.[2003). The uncer- 
tainty in the galactic B/D mass ratio is partly due to 
differences in the modelling of the disk component, where 
disk profiles exhibiting a central hole or depletion lead to 
B/D ratios at the high end, while double exponential pro- 
files without hole favour B/D ratios at the low end. Since 
we employed in our analysis disk models with central holes 
from Robin et al. ( 2003|) , we should for consistency com- 
pare our 511 keV B/D luminosity ratio to their (large) 
B/D ratio of ~ l. 6 But even with such large B/D mass 
ratios, the source population we are seeking for should still 
be at least 3 times more abundant in the bulge than in 
the disk of the Galaxy. We therefore conclude that the pri- 
mary positron source of the Galaxy is clearly associated 
with the galactic bulge. It therefore should belong to the 
old stellar population. 

Furthermore, the fact that the 511 keV emission 
matches well the morphology of the stellar bulge suggests 
that positron diffusion probably plays only a minor role. 
Were positron diffusion to be important we would expect 
to find substantial 511 keV emission in gas- rich regions 
adjacent to the rather gas-poor galactic bulge, such as the 
molecular ring structure at galactocentric distances of ~ 4 
kpc. However, we do not find any evidence for 511 keV 
emission correlated with this structure. We therefore con- 
clude that positron diffusion is negligible at galactic scales 
(i.e. kpc scales). 

5.3. Constraints on the disk source 

One certain source of positrons in the disk of the Galaxy is 
the radioisotope 26 A1. It decays with a lifetime of r ~ 10 6 
yr with emission of a 1809 keV gamma-ray photon; ~ 85% 
of the decays are also accompanied by the emission of a 
positron. The galactic distribution of 26 A1 is well known 

6 We note that, independently, large B/D ratios are also 
favoured by microlensing surveys towards the galactic bulge 
region (Binney & Evans 20011. 



thanks to observations of the COMPTEL telescope aboard 
CGRO, and follows that of the young stellar population. 
Thus, under the assumption that the positrons annihilate 
close to their production site, 511 keV line emission along 
the galactic plane is expected, showing the morphological 
characteristics of a young stellar population. 

The expected 511 keV line flux F 5n due to 26 A1 de- 
cay is related to the 1809 keV line flux Fi 80 g through 
F 5U = 0.85 x (2 - 1.5/p) x Fi 809 . Using the COMPTEL 
measurement of the 1809 keV flux along the galactic plane, 
Fi 809 = 9 x 10~ 4 ph cm^s" 1 (c.f. Table 4.3 in Knodlseder 
1997), and assuming / p = 0.93 leads to an expected 
511 keV line flux of F 5 n = 5 x 10~ 4 ph cm _2 s _1 . 

Fitting our model of the young stellar population 
(model DO), together with bulge models to the data sug- 
gests a disk flux in the range (4 — 6) x 10~ 4 ph cm _2 s — 
To explore the sensitivity of the disk flux on the assumed 
disk model we also paired the shell model with tracers of 
1809 keV line emission, such as the DMR free- free and the 
DIRBE 240 /im emission maps (Knodlseder et al. 1999b). 
This resulted in slightly larger (and more significant) disk 
fluxes of (8.3 ±2.3) x 10 -4 ph cm~ 2 s~ 1 . Comparing these 
values to what we expect from 26 A1 suggests that 60—100% 
of the galactic plane emission may be attributed to (3 + - 
decay of 26 A1. 

If this contribution is subtracted a 511 keV disk flux of 
at most ~ 3 x 10~ 4 ph cm~ 2 s _1 remains to be explained 
by other positron sources. A possible candidate is the ra- 
dioisotope 44 Ti whose daughter isotope 44 Sc decays via 
/3 + -decay into stable 44 Ca (r ~ 87 yr). In contrast to 26 A1 
there is no firm measurement of the present day galac- 
tic 44 Ti mass (nor of its spatial distribution), but simple 
chemical evolution arguments lead to the expectation that 
about 4 x 10~ 6 Mq of 44 Ti are produced per year (Leising 
& Share 1994). Under the assumption that all positrons 
escape the production site this yield translates into an 
annihilation rate of 3 x 10 42 s _1 . Assuming further that 
44 Ti is distributed following model DO, a 511 keV disk 
flux of ~ 8 x 10 -4 ph cm _2 s _1 is expected. In view of the 
approximate estimation and in particular in view of the 
uncertainty about the spatial distribution this value seems 
in reasonable agreement with the remaining 511 keV flux 
of - 3 x 10~ 4 ph cm^s" 1 . 

It is intriguing that the galactic disk flux could be en- 
tirely explained by the radioactive decay of 26 A1 and 44 Ti. 
This would suggest that once the 26 A1 and 44 Ti contri- 
butions have been subtracted only the bright bulge com- 
ponent of 511 keV emission remains, which would then 
demand a specific source population that is only con- 
fined to the inner Galaxy. However, we cannot immedi- 
ately draw this conclusion. Fitting an old stellar popula- 
tion disk model, with larger scale height and smaller radial 
scale than the young one, increases the estimate of the disk 
flux by about a factor of two, leaving room for a weak disk 
component not associated with 26 A1 or 44 Ti. This would 
suggest lower limits on the B/D flux (luminosity) ratio of 
*b 2 (i£ 6) for the source population that gives rise to the 
galactic bulge emission. 
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5.4. Constraints on the bulge source 

5.4.1. Massive stars, core collapse supernovae, pulsars 

The bulge dominance of the 511 keV emission immediately 
excludes scenarios in which the bulk of galactic positron 
production is related to massive stars. Such scenarios in- 
clude the production of /3 + -decay radioisotopes produced 
by Wolf-Rayet stars and all types of core collapse super- 
novae (including hypernovae) and the pair production in 
the strong magnetic fields of pulsars. Massive stars may 
well explain the faint disk component of 511 keV emission 
via the radioactive decay of 26 A1 and 44 Ti (c.f. section l5~3*|) . 
They cannot, however, explain the majority of the emis- 
sion, which would in that case resemble the 1809 keV line 
emission. 

5.4.2. Hypernovae 

Casse et al. (200 1 | proposed that a recent hypernova at 
the galactic centre could be responsible for the observed 
positron emission, but there is no observational evidence 
that such an explosion indeed took place. Hypernovae are 
believed to be related to Wolf-Rayet stars, which arc dis- 
tributed in the galactic disk following the spiral arm pat- 
tern, hence it would be much more likely to find a recent 
hypernova at an arbitrary position along the galactic plane 
(or at the tangent points of the spiral arms) rather than 
at the position of the galactic centre. And even in the 
rare event of a hypernova exploding right at the galactic 
centre, it would be difficult to explain why the resulting 
511 keV annihilation radiation (arising from the /3 + -decay 
of freshly synthesised 56 Co) should reflect the stellar mor- 
phology of the galactic bulge. We therefore conclude that 
it is unlikely that galactic centre hypernovae are the source 
of the bulge positrons. 

5.4.3. Cosmic-ray interactions 

Interactions of cosmic-ray particles with the ambient in- 
terstellar medium may produce positrons, primarily via 
the N + p — » 7r + — > e + reaction channel (N stands for nu- 
cleus) . The effect of the cosmic-ray interaction is best seen 
in the GeV gamma-ray domain, and has been comprehen- 
sively mapped by the EGRET satellite aboard CGRO. 
The EGRET all-sky map shows dominant emission from 
the galactic plane, which follows a linear combination of 
various gas and dust tracers in the galaxy. Thus cosmic-ray 
interactions should also lead to disk dominated 511 keV 
emission, which is at odds with our observations. 

5.4.4. X-ray binaries 

Positron production in X-ray binaries may occur either as 
a result of 77 pair production in the luminous compact 
region around the compact object or due to nuclear in- 
teractions that may form excited nuclei that subsequently 
decay through the emission of positrons. Galactic black 



holes and microquasars, where the positrons are ejected 
in a relativistic jet with Lorentz factors of a few (Dermer 
& Murphv |2()()1[) . are the two leading candidates. 

X-ray binaries are separated into two classes, depend- 
ing on whether the donor is a high-mass star (HMXB) 
or a low-mass star (LMXB). The two classes show clearly 
different spatial distributions (Grimm et al. 12002(1 . HMXB 
are associated with the young stellar population and are 
primarily found in the galactic disk. Consequently they 
can immediately be excluded as the source of the bulge 
positrons. LMXB, in contrast, are strongly concentrated 
towards the galactic bulge, and are more promising source 
candidates. Among the 150 LMXBs listed in the catalogue 
of Liu et al. (|2001f> . more than 50% are observed towards 
the galactic bulge. Correcting for completeness, Grimm et 
al. (2002) find a B/D ratio of ~ 0.9 and a vertical scale 
height of 410 pc for the LMXB distribution. 

Formally, the LMXB B/D ratio is considerably below 
the value required by our 511 keV data, yet the large verti- 
cal scale height of LMXB could lead to a scenario where a 
substantial fraction of positrons from disk LMXB may es- 
cape into the galactic halo. This scenario works as follows. 
Since the scale height of LMXB (410 pc) considerably ex- 
ceeds the scale height of the dense interstellar gas layer of 
our Galaxy (~ 100 pc), positrons from disk LMXB are 
ejected into rather low-density regions, typically a fac- 
tor of 10 — 100 less dense than regions found near the 
galactic plane fFerriere ll998[) . Before positrons can annihi- 
late they have to slow down considerably, mostly through 
Coulomb interactions, with a characteristic timescale of 
t"sd ~ lO 5 ?^ 1 yr, n being the ISM density in units of 
cm~ 3 (Forman et al. I198fi|) . Consequently, positrons live 
10 — 100 times longer at large scale heights than near 
the galactic plane, allowing for substantial diffusion be- 
fore annihilation takes place. The typical diffusion length 
depends much on the magnetic field configuration and the 
amount of ISM turbulence at large scale heights, but qual- 
itatively it seems plausible that a considerable fraction of 
the disk positrons may annihilate in the galactic halo. 

The resulting broad diffuse component of 511 keV 
emission would be difficult to detect with SPI. In par- 
ticular, the present dataset, for which good exposure is 
restricted to a band b ±15° along the galactic plane 
(c.f. Fig. [TJ, makes it virtually impossible to measure 
a disk component with a broad-latitude distribution of 
511 keV emission. Therefore it would be sufficient that 
~ 2/3 of the positrons produced by disk LMXB escape 
into the galactic halo to reconcile the LMXB distribution 
with a B/D ratio of 3. If we require a more extreme B/D 
ratio of 6, as expected after subtraction of the 26 A1 com- 
ponent from the disk (c.f. section f5.3|l . a positron escape 
fraction of ~ 80% would be needed for disk LMXB. It re- 
mains to be seen whether such large escape fractions are 
feasible. 

An alternative way to test the LMXB scenario is to 
search for 511 keV line emission from individual bright 
and/or nearby objects. So far no emission is seen towards 
the interesting candidates Sco X-l (the brightest LMXB) 
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and Cen X-4 and A0620-00 (probably the most nearby 
LMXB at ~ 1.2 kpc), but we plan for deep observations 
of these objects in the near future to search for their anni- 
hilation signatures. The detection of positron annihilation 
signatures from nearby objects could however be ham- 
pered by (even modest) positron diffusion away from the 
sources, which would lead to extended 511 keV emission 
halos around the objects. So even for the modest angular 
resolution of SPI of ~ 3° nearby individual LMXB could 
appear as extended sources, and their low surface bright- 
ness could make their detection more difficult. 

5.4.5. Classical novae 

Among all proposed positron candidate sources, classi- 
cal novae, i.e. thermonuclear runaways on white dwarfs 
in accreting binary systems, are the sources for which 
the largest B/D ratios of ~ 3 — 4 have been suggested 
(Delia Valle & Duerbeck 113951 Delia Valle & Livio ll()94a|) . 
Interstellar extinction, in particular towards the galactic 
bulge region, makes it virtually impossible to derive their 
spatial distribution in the Milky Way directly, but no- 
vae are readily observed in nearby external galaxies which 
may serve as templates (e.g. Shatter et al. .2000 ) . Due to 
its proximity and due to its similarity to the Milky Way, 
M31 is the primary source of information, and modern 
investigations indicate that novae reside primarily in the 
bulge region of M31 (Ciardullo et al. 119871 Capaccioli et 
al.riTIs!)). Although it had been suggested that selection ef- 
fects may have "faked" such a finding (Hatano et al.[l997) 
the recent study of Shatter & Irby l|2001[l demonstrates 
that such biases, if they exist, must be small. 

Novae produce positrons via the /3 + -decay of radioac- 
tive isotopes synthesised during the thermonuclear run- 
away, mainly of 13 N, 18 F, and 22 Na (lifetimes r = 14 min, 
2.6 hr, and 3.75 yr, respectively). 22 Na yields of 6 x 10~ 9 
Mq, as suggested by theoretical nucleosynthesis calcula- 
tions for ONe novae (Hernanz et al. I2002f) . would require 
nova rates of ~ 1600 yr -1 to maintain positron production 
and annihilation in an equilibrium state, a rate which is 
considerably above estimates of 35 ± 11 yr" 1 for all types 
of novae in the entire Galaxy ( Shatter HM7| . 13 N yields 
of 2 x 10 -7 Mq for low-mass CO novae are more promis- 
ing (Hernanz et al. 12002) since they would require nova 
rates of only 26 yr" 1 if all positrons could indeed escape 
from the nova envelope into the ISM. However, with a 
13 N lifetime of 14 min it seems unlikely that this would 
be possible. 

It is probable that large fractions of the 13 N positrons 
annihilate within the dense nova envelope, leading to 
prompt annihilation that could give rise to transient an- 
nihilation signatures (Leising & Clayton 119871 Gomez- 
Gomar et al. I1998fl . This signature has been sought using 
various gamma-ray telescopes, but has so far eluded de- 
tection (see Hernanz & Jose 2004 and references therein). 
Detection of the transient signature may help to shed light 
on the positron escape fraction, and could show whether 



novae contribute to the galactic bulge positron budget or 
not. 

5.4.6. Thermonuclear supernovae 

In view of their potential to produce large numbers of 
positrons, thermonuclear Type la supernovae (SN la) are 
often considered as the most plausible source of positrons 
in the Milky Way (Dermer & Murphy IM)T|l . SN la pro- 
duce positrons via the /3 + -decay of radioactive 56 Co (r = 
111 days). Expected 56 Co yields of ~ 0.6 Mq provide 
~ 2.5 x 10 54 positrons per event, although, as with novae, 
prompt annihilation in the supernova envelope probably 
prevents large fractions of the positrons from escaping into 
the ISM. From the analysis of late light curves of SN la 
Milne et al. I(1999j) derive a mean escaped positron yield 
of ~ 8 x 10 52 positrons per SN la, corresponding to a 
positron escape fraction of / ~ 0.03. A recent study of 
SN 2000cx even suggests / ~ 0, but SN 2000cx was an 
unusual event that may not represent the average SN in 
the bulge of our Galaxy (Sollerman et al. 2004). 

Assuming therefore / = 0.03 a bulge SN la rate of 0.6 
per century is required to maintain the observed 511 keV 
luminosity in a steady state. Unfortunately the rate and 
distribution of SN la in our Galaxy are only poorly known. 
The galactic SN la rate is generally inferred from rates 
observed in external galaxies which are then scaled to the 
mass and the type of the Milky Way. In that way rates 
of 0.3 — 1.1 SN la per century are derived (Tammann et 
al. ESI| Cappellaro et al. HM7I Mannucci et al. EHEl, 
sufficient to maintain the galactic 511 keV luminosity. In 
contrast, when we follow the suggestion of Prantzos ( 2004 ) 
and derive the bulge SN la rate by scaling the SN la rate 
observed in early-type galaxies to the mass of the galac- 
tic bulge, a bulge SN la rate of 0.08 SN la per century is 
obtained. This value is much too low to explain the ob- 
served bulge 511 keV luminosity. It is difficult to judge if 
the galactic bulge can indeed be considered as a downsized 
version of an elliptical galaxy, in particular in view of the 
differences in the evolution of the galactic bulge and an 
elliptical galaxy. Furthermore, it is suggested that differ- 
ent SN la explosion mechanisms exist in different types 
of galaxy (e.g. Delia Valle & Livio I1994bl Howell IMTTl 
Mannucci et al. 2005 ) making the proposed extrapolation 
even more uncertain. 

Observations of external galaxies indicate that SN la 
distributions are strongly peaked towards galactic centres 
(Bartunov et al. H992fl . yet reliable determinations of B/D 
ratios are difficult in view of observational biases and selec- 
tion effects (Wang et al. IT537I Hatano et al. ITMSI Howell 
et al. I2000|l . If the SN la distributions follows that of no- 
vae (both populations are believed to arise from accreting 
white dwarfs; see van den Bergh ll988fl one can expect B/D 
ratios of ~ 3 — 4. Even higher B/D ratios can be achieved if 
part of the positrons produced by disk SN la, which have 
a vertical scale height of ~ 330 pc (Chen et al. 1200 ill, es- 
cape into the halo (c.f. section IB. 4. 4(1 . Thus SN la could 
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indeed present the required characteristics that are needed 
to explain the positron distribution and annihilation rate 
in the Galaxy. 

Alternatively, instead of explaining the bulge emis- 
sion globally we may search for 511 keV emission from 
nearby Type la supernovae remnants, such as SN 1006, the 
Lupus Loop, or the Tycho SNR. Assuming a mean escaped 
positron yield of ~ 8 x 10 52 positrons per SN la (Milne 
et al. I1999f) and a positronium fraction of / p = 0.93, the 
mean expected 511 keV line flux from an individual SN la 
is estimated to 

F 5U = 1.3 x ID"* (^p) 2 (^) P h cm-s- (9) 

where D is the distance to the supernova remnant in units 
of kpc, and r is the mean lifetime of the positrons in years, 
between 10 3 — 10 7 yr, depending on the density, temper- 
ature, and ionisation state of the annihilating medium 
(Guessoum et al. 119910 . Taking distances of D ~ 2 kpc 
for SN 1006 (Laming et al. HMfy . of D ~ 1 kpc for the 
Lupus Loop (Leahy et al. I1991[) . and of D ~ 2.3 kpc 
for Tycho (Hughes 2000) and assuming a positron life- 
time of 10 5 yr results in predicted 511 keV line fluxes 
of 0.3 x 10~ 4 ph cm^s" 1 , 1.3 x 10~ 4 ph cm^s" 1 , and 
0.2 x 10 -4 ph cm~ 2 s _1 , respectively. At least for the Lupus 
Loop, the predicted flux is close to our upper 511 keV 
flux limit, indicating that dedicated deep observations of 
nearby supernova remnants can help to answer the ques- 
tion about the galactic positron source. Such dedicated 
observations with INTEGRAL are already scheduled. 

5.4.7. Light dark matter annihilation 

Light dark matter (1-100 MeV) annihilation, as suggested 
recently by Boehm et al. (2004), is probably the most ex- 
otic but also the most exciting candidate source of galactic 
positrons. Unfortunately, the spatial distribution of dark 
matter in general, and light dark matter in particular, is 
only poorly constrained by observational data, at least for 
the inner Galaxy. The debate of whether the dark mat- 
ter profile shows a cusp towards the galactic centre is still 
not settled, but it seems clear now that, dynamically, dark 
matter plays only a minor role in the inner 3 kpc of the 
Galaxy. In this region the stellar mass dominates (Binney 
& Evans IM7T1 Klypin et al. TO^ . 

Due to these uncertainties it is difficult to judge 
whether the observed 511 keV emission could be explained 
by dark matter annihilation. Maybe more promising is the 
idea to search for signatures of dark matter annihilation in 
nearby, external galaxies. Hooper et al. I|2004|l suggested 
that nearby dwarf spheroidal galaxies may provide promi- 
nent sources of 511 keV line emission due to the high 
densities of dark matter that are known to be present. 
They proposed the nearby Sagittarius dwarf galaxy (Sgr 
dwarf) as most promising candidate and estimate 511 keV 
line fluxes in the range of (1 — 7) x 10~ 4 ph cm _2 s _1 , de- 
pending on the assumed dark matter halo profile. Cordier 



et al. H2004D searched for emission from this Galaxy us- 
ing SPI and obtained a (3cr) upper limit of 3.8 x 10~ 4 
ph cm _2 s _1 . Our upper limit of 1.7 x 10~ 4 ph cm^s" 1 is 
substantially lower, and excludes almost all types of halo 
models for this galaxy, in particular those with a central 
cusp. 

Standard cold dark matter cosmology predicts cuspy 
dark matter distributions (Klypin et al. I2002|l . so in prin- 
ciple Sgr dwarf should have been detected by SPI if dark 
matter annihilation were a viable scenario. Maybe dark 
matter halos are less cuspy than theory predicts? This 
possibility is indeed indicated by observations of our own 
Galaxy (Binney & Evans 2001) and dwarf galaxies (Blais- 
Ouellette et al. 119991 Kleyna et al. 1277113")) . But in this case 
dark matter annihilation should not lead to a compact but 
to a rather extended 511 keV emission feature - in con- 
tradiction to what SPI observations of the inner Galaxy 
suggest. From the arguments given one may question the 
dark matter scenario. However, it is certainly premature 
to reject them totally because of the uncertainties in the 
dark halo profiles and the annihilation conditions. 

6. Conclusions 

Our first mapping of 511 keV gamma-ray line emission 
over a large fraction of the celestial sphere leads us to the 
following observations: 

1. 511 keV emission is significantly (~ 50cr) detected to- 
wards the galactic bulge region, and, at a very low level 
(~ 4<r), from the galactic disk 

2. there is no evidence for a point- like source in addition 
to the diffuse emission, down to a typical flux limit of 
~ 10~ 4 ph cm _2 s _1 

3. there is no evidence for the positive latitude enhance- 
ment that has been reported from OSSE measure- 
ments; the 3<T upper flux limit for this feature is 
1.5 x 10~ 4 ph cm~ 2 s _1 

4. the bulge emission is spherically symmetric and is cen- 
tred on the galactic centre with an extension of ~ 8° 
(FWHM); it is equally well described by models that 
represent the stellar bulge or the halo populations 

5. the bulge annihilation rate is (1.5± 0.1) x 10 43 s^ 1 , the 
disk annihilation rate is (0.3 ± 0.2) x 10 43 s _1 

6. the bulge-to-disk luminosity ratio lies in the range 3 — 9 

The bulge dominated 511 keV line emission morphol- 
ogy suggests an old stellar population as the main galactic 
positron source. In contrast, the faint disk emission is well 
explained by the release of positrons during the radioac- 
tive decay of 26 A1 that originated from massive stars, with 
a possible contribution from 44 Ti synthesised during su- 
pernova explosions. 

The extreme bulge-to-disk ratio that is observed in 
the 511 keV luminosity imposes severe constraints on the 
principal galactic positron source. Type la supernovae, 
low-mass X-ray binaries or dark matter annihilation may 
possibly satisfy these constraints, but uncertainties in the 
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knowledge about the spatial distribution of these objects 
and the positron escape processes prevents us from draw- 
ing firm conclusions. Novae could probably most easily 
explain the large B/D ratios, yet an implausibly large 
positron escape fraction from 13 N decay would be re- 
quired to accommodate the observed annihilation rate. 
SN la could explain the annihilation rate for a modest 
positron escape fraction, but it is questionable if they have 
the required large B/D ratio. LMXB could reproduce the 
observed B/D ratio provided that a substantial fraction 
of positrons ejected by disk LMXB escape into the halo. 
Light dark matter is an exciting option, but it remains 
to be seen if the observed 511 keV emission distribution 
is compatible with the profile of the galactic dark matter 
halo. 

Future deep observations of individual nearby candi- 
date sources may provide the means to identify the galac- 
tic positron source. As such we will soon observe the X-ray 
binaries Sco X-l and Cen X-4 with INTEGRAL, and ob- 
servations of the nearby supernova remnant SN 1006 are 
already scheduled. We cannot be sure that any of these ob- 
servations will allow the detection of a 511 keV signal, but 
were such a signal detected we would gain important new 
insights in the primary source of positrons in our Galaxy. 

Appendix A: Summary of galactic density profiles 

Bulge models: We model the galactic bulge as a triax- 
ial stellar bar for which the parameters are summarised 
in Table IA.1I The apparent intensity distribution on the 
sky is computed in the galactic frame which we define 
by a right-handed cartesian coordinate system where the 
Sun is located on the negative y-axis, at —Rq = 8.5 kpc. 
The transformation from the galactic frame into the bar 
frame is performed by two consecutive rotations: the first, 
represented by the matrix Rx(ct), is a counterclockwise 
rotation by an angle a around the z-axis; the second, rep- 
resented by the matrix i?y (—/?), is a clockwise rotation by 
an angle j3 around the new y-axis, i.e. 



r' = R Y (-/3)Rx(a)r 

The effective bar radius is defined by 



i? s = 




where 

r i.o 

is a cutoff function and 



(A.l) 



(A.2) 



where a x , a y , and a z are the scale lengths and C± and C|| 
are the face-on and the edge-on shape parameters. The 
radial dependencies of the bar density are given by various 
density profiles that we designate by labels: 

(iky) x exp(-0.5i? s 2 ) 



G 
G3 

E 
E3 

S 

P 



Pg = Po x /„ 

PO X fn 
PO X fn 
PO X fn 
PO X fn 



PG3 
PE -- 

PE3 
PS - 



c(-Rxy) 
c(-Rxy) 
c(-Rxy) 
c(-Rxy) 



xi?- 18 exp( 
x exp(-iC) 
x K {R S ) 
x sech 2 (i? s ) 



PP = PO X /max(i?xy) X [1 + (R s / R c )\ 



(A.3) 
(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 



for R XY < i? ma>: 



flma= 



for R xy > Rmax 



i?xy - V* 2 + y 2 (A. 10) 

is the distance from the galactic centre in the xy-plane. 

Model DO: For the young disk population we use the 
model proposed by Robin et al. H2003[) to describe the 
young (age < 0.15 Gyr) stellar disk population of the 
Galaxy: 

p(x, y, z) = p (cxp(-(a/i? ) 2 ) - exp(-(a/i? 2 ) 2 )) (A.ll) 
where 



2 i 'I i 111 

x + y + z /e 



(A.12) 



This model presents a truncated exponential disk profile 
with a fixed disk axis ratio of e = 0.014, a fixed disk scale 
radius of Ro — 5 kpc, and a fixed inner disk truncation 
radius of Ri = 3 kpc. The vertical exponential scale height 
of the disk is zq = 70 pc. 

Model Dl: For the old disk population we use the model 
proposed by Robin et al. (2005J) to describe the old (age 
7—10 Gyr) stellar disk population of the Galaxy: 

p(x,y,z) = p Q (cM-(0-^ + a 2 /R 2 Q ) 1/2 )- 

exp(-(0.25 + a 2 /i? 2 ) 1 / 2 )) (A.13) 

This model presents a truncated exponential disk profile 
with a fixed disk axis ratio of e = 0.0791, a fixed disk scale 
radius of Rq — 2.53 kpc, and a fixed inner disk truncation 
radius of Ri — 1.32 kpc. The vertical exponential scale 
height of the disk is z$ = 200 pc. 

Model H: To model the stellar halo we use the general 
model proposed by Robin et al. (|2003|l : 



p(R,z)=p (a/R Q )- n (A.14) 

where a is defined by Eq. IA.12I (the flatness of the model 
is set by the value of the axis ratio e in Eq. IA.12f> . and 
a > a c , avoiding a singularity at the galactic centre, n 
determines the slope of the density profile. 

Model S: A set of galactocentric nested spherical shells of 
constant density. The radii of the shells and their number 
has been varied in order to maximise the MLR, while using 
the minimum number of shells required to satisfactorily 
describe the data. 
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Table A.l. Parameters of the triaxial galactic bulge models used in this work (see text). 



Model 


a 


P 


Ox 


a y 


a z 


c± 


C\\ 


n 


Rc 


Rmax 


Q>max 


GO 


0.0 


0.0 


0.91 


0.51 


2 


2 






OG 




Gl 


73.8 


0.6 


3.11 


0.76 


0.48 


2 


2 






2.40 


0.71 


G2 


70.7 


0.8 


1.56 


0.60 


0.45 


2 


4 






2.40 


0.71 


G3 


59.0 


0.7 


3.78 


1.44 


1.19 


2 


2 






2.40 


0.71 


E0 


0.0 


0.0 








2 


1 






oc 




El 


65.7 


0.5 


2.23 


0.65 


0.32 


1 


1 






2.40 


0.71 


E2 


49.0 


0.6 


0.75 


0.19 


0.27 


2 


2 






2.40 


0.71 


E3 


49.8 


0.7 


0.69 


0.19 


0.28 


2 


4 






2.40 


0.71 


Sf 


76.2 


0.02 


1.70 


0.64 


0.44 


1.57 


3.50 






3.13 


0.46 


Ef 


80.5 


-0.05 


1.89 


0.66 


0.43 


1.61 


3.49 


1.44 




3.57 


0.56 


Pf 


76.8 


0.02 


1.81 


0.65 


0.43 


1.65 


3.02 


5.04 


1.23 


2.71 


0.88 


SpR 


79.4 


-0.8 


1.87 


0.56 


0.47 


3.49 


3.38 






3.71 
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